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The quantum mechanical probability densities are compared with the probability densities treated 
by the theory of random variables. The relevance of their difference for the interpretation of quantum 
mechanics is commented. 
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I. INTRODUCTION 

Probabilities were introduced in quantum mechanics by Max BornjlJ when he proposed an interpretation for the 
wave function ^(x) of a particle as a probability amplitude. Thereby ^(a;)) 2 is a probability density, that is, the 
probability assigned to the event corresponding to the location of the particle in the interval [x,x + dx]. Today, 80 
years later, we don't know the nature of this probability, that is, we don't know whether this probability has an 
ontological or a gnoseological character: we don't know whether the particle is located somewhere but we can not 
. know where it is, and the best that we can do is to give a probability for it, or on the contrary, the very location of the 
particle is diffuse by nature. Although the Bell-Kochen-Specker theoremP,Q favours an ontological interpretation, 
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as can be easily seen in a recent illustration [J] of the theorem concerning position and momentum observables, the 
opposite view, assuming context dependent definite values for the observables, is not excluded. It soon became clear 
however that these probabilities do not behave as such for mutually exclusive events and the rule "thou shall not add 
the probabilities but their amplitudes instead" was adopted although never understood. A beautiful illustration of 
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, the necessity of this rule is given by R. Feynman|5J in the analysis of the two slits experiment. 

Quantum mechanical "probabilities" are not added; however, probabilities are studied by a rigourous mathematical 
theory of random variables that tell us that, indeed, they can be added. Instead of requiring a change in the 
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mathematical theory we should perhaps accept that the quantum mechanical densities are not probabilities. In this 
work we will analyse the quantum mechanical densities and compare them with the probability densities defined in 
the theory of random variables. We will see that, in general, the quantum mechanical densities do not behave as the 



probability densities of random variables and we will comment on the consequences of this for the interpretation of 
quantum mechanics. 
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II. RANDOM VARIABLES 



For completeness, we will recall some concepts of the theory of random variables. A random variable A is a 
mathematical entity that may be assigned numerical values a with some probability. We can use the symbol A => a 
to denote the event that the random variable is assigned a value in the interval [a, a + da]. We define the Probability 
Density Function g(a) as the probability of realization or appearance of the event A => a. Of course, theses densities 
are nonnegative and normalized. 

Probabilities play an important role in classical physics and are omnipresent in quantum mechanics so, in some cases, 
it is interesting to consider the obscrvablcs of a system as random variables with their corresponding probability density 
functions. For instance, we may study the possibility of assuming that the position of a particle X and its momentum 
P are random variables with their associated probability densities g(x) and 7r(p) that give the probability of location 
of the particle and of having some value of momentum. Since position and momentum are two independent concepts, 
in the sense that a particle at a given position can have any value of momentum, and some value of momentum can be 
realized at any position, we can assume that position and momentum are described by independent random variables. 
We will later mention the possibility that they are not independent. 

Given a random variable A with density g(a), a function H (A) of the random variable A, will have a density rj(h) 
given by 



This result is intuitive: the probability for some value h is the probability of the event A =>■ a times the probability that 
h = H(a), which is a sharp Dirac distribution, added for al possible values of a. Similarly, it is easy to understand that 
if A and B are independent random variables and g(a) and ir(b) are the corresponding probability densities associated 
with the events A => a and B => b, then the probability density cr(c) assigned to the random variable C = F(A, B) is 
given by 



This result can not be avoided. If by any reason, the event C => c doesn't have the given probability density, it does 
not means that the theory is wrong but it suggests that, perhaps, g(a) and ir(b) are not the probability densities of 
the random variables A and B. 



We will prove the following theorem. Assume two observable of a quantum system represented by the operators 
A and B in a Hilbert space. Let g(a) and w(b) be the distributions of the eigenvalues of the operators as predicted 
by quantum mechanics for the system in a state that is, g(a) = (^>,P a ^/) where P a is the projector in the 
invariant subspace of the eigenvalue a (in the case of nondegeneracy of the eigenvalue a with eigenvector ip a , we 
have P a — (p a ((p ai u) and we get the more familiar expression g{a) = \(<p a , ^)| 2 ) an d similarly, 7r(6) = (^>,Pb^>) 
(n(b) — \ ((j>b, ^)\ 2 if b is nondegenerate with eigenvector Let C — F(A, B) be an observable built as a function of 
the observables A and B whose operator has the eigenvalues c and the corresponding projector P c (eigenvector £ c if 
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III. QUANTUM RANDOMNESS 



nondegenerate). Then, considering the quantum mechanical distribution ctqm(c) and the probability density provided 
by the theory of random variables <jrv(c) given by 



a QM (c) = <*,P C *> (3) 
<j RV {c) = jdajdb g(a) n(b) S(c- F(a,b)) , (4) 

we can prove the following: 

1. if [A, B] ^ then in general it is ctqm(c) ^ ctrv(c); 

2. if [A, B] = then we have two cases: 

(a) if they commute because they are related by a function like B = G(A), then ctqm(c) = <trv(c)', 

(b) if they commute because they are independent observables related to different degrees of freedom, then we 
have again two cases: 

i. if the state is factorizable in the product Hilbert space corresponding to the different degrees of freedom, 
then ctqm(c) = ctrv(c); 

ii. if the state is not factorizable (entangled) in the product Hilbert space corresponding to the different 
degrees of freedom, then ctqm(c) ^ <jrv{c). 



Summarizing we have 



[A,B}^0, ctqm{c)^<j rv {c); 

B = G(A), <tqm(c) = <7 RV (c) 



[A,B]=0, { 



A <g)l,l(8)S, 



(5) 

* = x ® "0, ctqm(c) = <7rv(c); 

* X ® "0, o-qm(c) 7^ cri?y(c)- 



That is, in general the quantum mechanical distribution and the random variables densities are different and only in 
some cases (commuting operators and factorizable state) they can be equal. 

In order to compare them, we can write Eqs.(3) and (4) in a way to emphasize their difference (assuming nonde- 
generate eigenvalues a, b and c) as 

vqm(c) = Jdajdb (<fo, &(&,•) <Pa) (6) 

<Jrv{c) = [dafdb <y(c-F(a, &))*<*,•) y> ) , (7) 



The QM case involves the projector in the state whereas in the RV case it involves the projector in the state 
and a constraint S (c — F(a, b)). Of course this formal difference only suggests that (Jqm ^= vrv but it docs not prove 
it because the different operators are inside an integral. We will later prove anyway that they are indeed different. A 
nice comparisons of the two distributions (3) and (4) is obtained when we calculate the moments of the distributions. 
For the QM case we have 



E ( qI = Jdc c n <j QM (c) = (*,Jdc c n P c ^ = (tt.C"*) = ($,r(A5)$> , 



(8) 



and in the RV case we have 

E$ = Jdc c n o RV (c) = Jdc c n Jdajdb g{a) n(b) S(c- F(a, b)) = Jdajdb F n (a, b)g(a) %{b) . (9) 

Inserting the completeness relation in Eq.® and replacing the quantum mechanic expression for the densities g and 
7r we get 

E QM = J da J db (*,PaF n (A,B)P b V) , (10) 

E RV = J da J db (y,PaF n (a,b)P*P b y) . (11) 

Notice that the QM case involves the operator F n (A,B) and the RV case the number F n (a,b) multiplied by the 
projector in the state \I>. A nicer comparison is obtained in the particular case where [A, B] = 0. In this case we can 
write F n (A, B) as a power series with the operator B always to the right and A to the left and using the relations 
BP b = bP b and P a A = P a a we have 

u - =n _ . 4m = Jdajdb F n (a,b)(V, P a P b *) = Jda Jdb F n (a,b) lr{P a P b P,} , 

^ E$ = Jda Jdb F n {a,b) (^,P a P, v P b *) = Jda Jdb F n (a,b) Tr {P a P^P b P*} , 

where we have used the trace expression for the expectation values. Clearly the moments, and therefore the distribu- 
tions, are different in general although for some states they may be equal (for instance when P b and P^ commute). 

Considering some examples it is clear that for noncommuting observables, like position and momentum, it must be 
o~qm 7^ o-RV- Take for instance A = X 2 , B = P 2 and C = X 2 + P 2 where g(a) and 7r(6) are continuous, and therefore 
&rv(c) is a l so continuous, whereas o~qm{c) is not continuous because it corresponds to the energy distribution of an 
harmonic oscillator. 

Another case of physical interest with noncommuting observables is when A = X, the position of a free particle 
at t = 0, B = ^P, its velocity multiplied by the time t, and C = A + B = X + P-^ corresponds to the position of 
the particle at time t. The quantum mechanical prediction for the distribution of position at time t clearly differs 
from the random variable distribution. To see this we can, for instance, calculate the width of both distributions 
A 2 = E^ - (E^) 2 using Eqs.(8,9), or directly from Eqs. (3,4), and we obtain 

A 2 QM = A 2 x + ^A 2 p + -((y,(XP + PX)y)-2(y,Xy){y,Py)) , (13) 
A 2 RV = Al + — 2 A\ , (14) 

where A^ and A p are the width of the probability densities of position and momentum g(x) and 7r(p). The conceptual 
relevance of this difference will be discussed later. Notice that even if the observables A and B would commute, we 
would still get different results in some states where the corresponding correlation term in Eq.(13) does not vanish. 

Let us consider now the case of [A, B] = because B — G(A). In this case, if we know the density g(a) for the 
random variable A then we calculate the density of B using Eq.{T]) as 



tt(6) = Jda g(a) 5 {b - G(a)) . (15) 

This result of the theory of random variables is also obtained in quantum mechanics where n(b) = {^,P b ^). If 
[A, B] — 0, the operators A and B share eigenvectors, that is, they have the same projectors. However the function 



G(A) can introduce more degeneracy because there may be different values a ^ a' with G(a) = G(a') corresponding 
to the same eigenvalue b (think for instance that for X 2 we have P x i = P x + P-x)- Therefore the projector P b is 
obtained adding all projectors P a where a and b satisfy the condition b — G(a). That is, 

P b = Jda P a S{b- G(o)) . (16) 

Now, taking the expectation value P,*} we immediately obtain the relation in Eq. (fT5)) . In this case, the observable 
C = F(A,B) does not really depends on A and B but is a function of just one observable C — F(A,G(A)) and 
repeating the same argument it is clear that the theory of random variables and quantum mechanics predict the same 
density: (*, P c #) = Jda g(a) 5 (c - F(a,G(a))). 

We can now analyse the case where the observables A and B commute because they correspond to different degrees 
of freedom. These operators act in different Hilbert spaces H\ and H2 and the system is described by a state in the 
space H = Hi ® H-2- Given a state 4", there is a basis {xk} in Hi and a basis {C/c} in H2 that allow the bi-orthogonal 
decomposition of the state 

N 

* = ak * k ® ^ ■ 

k=l 

As is well known, if N = 1 the state is factorizable and if N > 2 the state is entangled or nonfactorizable. The operators 
and their corresponding projectors are extended in the product space as A® 1, P a ® 1, 1 ® B, 1 ® Pi,, and the densities 
associated with the observables are given by quantum mechanics as g(a) — ("P, P a <8> and 7r(6) = (\&, 1 <£> PfcW). 
The operator C = F(A, B) will have projectors P c associated with the eigenvalues c that project in a subspace of 
Hi (8 Wa- We will prove that 



P c = Jdajdb 6(c- F(a,b)) P a ®P b . (18) 
In order to prove it, let us consider the function F(A, B) expanded as a formal power series 

C = F(A, B) = J2 Cn, m A n ® Bm ■ (19) 

n,m 

Now, from the spectral decomposition of A = Jda a P a and using the orthogonality and idempotent property of the 
projectors P a P a ' = S( a ~ a ')Pa we easily find that A n = Jda a n P a and similarly B m = Jdb b m Pi,. Replacing above 
and reconstructing the function we get 



C = F(A, B) = jdajdb F{a, b) P a ®P b . (20) 

Now we write the function F(a, b) as an integral over c with a Dirac distribution 

C = F(A, B) = Jdajdb JdccS(c-F(a,b)) P a ®P b = Jdcc Jdajdb 6 (c - F '(a, b)) P a ® P b . (21) 

But this is precisely the spectral decomposition of the operator C and therefore the double integral over a and b is the 
projector P c as given in Eq. (|18|1 . We obtain the quantum mechanical distribution of the eigenvalues of the operator 
C by taking the expectation value in Eq. lfTH]) 

vqm(c) = (% P c *) = fdafdb S (c - F(a, b)) (% P a ® P b *) . (22) 



In this expression we consider now the two cases for the state, factorizable or entangled. In the first case the state is 
given by = \ ® £ and inserting above we have 

<?qm(c)= Jdajdb 5 (c-F(a,b)) ( X , P a x)((, PbQ , (23) 

but (x, PaX) = (*, Pa ® 1*) = g(a) and «, P b () = ir(b) and therefore 

a QM (c) = jdajdb S(c- F(a, b)) g(a) tt(6) , (24) 

which is the density predicted by the theory of random variables; that is, ctqm(c) — &rv{ c )- In the non- factorizable 
case we have 



ctqm(c) = Jdajdb 5 (c- F(a,b)) y^a%a r {Xk,PaXr)((k,Pb(r) , (25) 

k.r 

which contains all non-diagonal terms. Notice that here g{a) — J2k \ a k\ 2 {Xk, PaXk) and 7r(6) = ^ fc |afc| 2 (0t, Pb(k) 
and then 

<trv(c) = Jdajdb 6 (c - F(a,b)) J2\ a k\ 2 \a r \ 2 (Xk,PaXk)(Cr,PbCr) , (26) 

k,r 

therefore in general we have <tqm{c) ^ o\rv(c). 



IV. CONCLUSION 

The general conclusion that we can draw from this study is that quantum mechanics is not a random process in 
space. If it where, then the randomness in quantum mechanics should be described by the theory of random variables 
and we have seen that this is not the case. 

In particular, we saw in the last section that if we know the probability density for the position of a particle g(x) 
at one time, and we know the probability density of its velocity (or momentum) w(p), then, according to quantum 
mechanics, the probability density for the position at a later time is not calculated by sampling one position and adding 
a displacement sampled from the velocity distribution, as would be calculated with the theory of random variables. 
Notice that this result contradicts the gnoseological interpretation of the probability for position and velocity of the 
particle. Indeed, if the particle has a definite position, unknown to us, with its value distributed according to g(x), 
and the particle has a definite velocity, unknown to us and distritbuted according to n(p), then the position at a 
later time must given by the rule "xo + vt" where xo and v are taken from their distributions. In other words, the 
location at time t must be given according to the theory of random variables; but this is contrary to the quantum 
mechanical result! This can only mean that X and P are not random variables and g(x) and w(p) are not densities of 
probabilities although they are, in principle, measured as if they were, that is, by counting the frequency of realization 
of the events X =4- x and P =4- p. To call "probability" to something that is not a probability, that is, something 
that does not complies with the rules of probabilities, is perhaps a misuse of language and it could be convenient to 
denote the quantum mechanical distributions with another name. It is unfortunate that the term "probability" is 
irreversibly installed in quantum mechanics because, strictly speaking, they are not probabilities and the misnomer 
introduces confusion, not only in the teaching, but also in search of an interpretation of quantum mechanics. A more 



appropriated name for g(x) and 7r(p) could be for instance the existential weight of the events X =>- x and P => p 
in the system in a particular state. Although these existential weights are measured in an experiment in the same 
way that one measures a probability, these quantities are not probabilities and do not necessarily obey all the rules of 
probabilities dictated by the theory of random variables. Probabilities propagate according to the theory of random 
variables but the existential weights are calculated by the rules of quantum mechanics. Accordingly we can think 
about the position of a particle as a diffuse observable without a putative value. Every possible value, every event 
X =>• x, has a propensity to appear in an observation given by the existential weight, but this does not means that the 
event has that value with some probability and that the measurement of the observable exhibits a pre-existent value. 
This idea has been synthetically expressed by A. Peres saying that "unperformed experiments have no result" 

The results presented in this work suggest that we can not consider g{x) and Tv(p) to be the probability densities 
of two independent random variables associated with the position and momentum of a particle. We should consider 
now the possibility that position and momentum are not independent and are related by some unknown connection 
that is responsible for the quantum mechanical correlations. We could then recover the theory of random variables for 
quantum mechanics if we could show that g{x) and 7r(p) are the marginal distributions of a joint probability density 
W(x,p). There are many possibilities to achieve this, the most famous is the Wigner phase space distribution, but 
it was shownQ that the most general function producing the same expectation values as quantum mechanics is not 
nonnegative everywhere and therefore it can not be considered as a genuine joint probability density. 
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